Coupling of excitons and defect states in boron-nitride nanostructures 
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The signature of defects in the optical spectra of hexagonal boron nitride (BN) is investigated 
using many body perturbation theory. A single BN-sheet serves as a model for different layered BN- 
nanostructures and crystals. In the sheet we embed prototypical defects such as a substitutional 
impurity, isolated Boron and Nitrogen vacancies, and the di-vacancy. Transitions between the deep 
defect levels and extended states produce characteristic excitation bands that should be responsible 
for the emission band around 4 eV, observed in luminescence experiments. In addition, defect 
bound excitons occur that are consistently treated in our ab initio approach along with the "free" 
exciton. For defects in strong concentration, the co-existence of both bound and free excitons adds 
sub-structure to the main exciton peak and provides an explanation for the corresponding feature 
in cathodo and photo-luminescence spectra. 

PACS numbers: 71.35.-y 71.55.-I 
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I. INTRODUCTION 

Boron nitride crystals and nano structures are excel- 
lent wide band gap materials for the realization of UV- 
lasersi 1 - Defects formed during the BN-synthcsis or intro- 
duced as impurities can bind excitons and act as centers 
for their recombination. This leads to (desired or unde- 
sired) luminescence bands and losses in the photon fre- 
quency range of interest. In fact, recent experiments 2 - - — 
have given evidence that defects strongly effect the lu- 
minescence of BN-based materials. In order to achieve 
optimal optical properties, the relevant defects have to 
be identified such that methods can be devised to tune 
their concentration in the material. 

Boron nitride (BN) is known to exist in three different 
crystals structures, namely hexagonal BN (h-BN), cu- 
bic BN (c-BN), and wurtzite BN (w-BN). Among these 
structures, h-BN is the stable one at room tempera- 
ture and ambient pressure. Like graphite, it consists 
of stacked BN-layers with a honeycomb structure. BN 
nanostructure o 10 ! 11 ranging from single and multi-wall 
nanotubes to BN-fullerens and nanocones can be envis- 
aged as having been tailored from a single BN sheet. 
Unlike in graphite, the partially ionic character of the 
BN bond results in a wide band-gap of about 6.5 eV for 
bulk h-BNi£. The combination of a large gap with a 
strong electron-hole attraction makes the optical proper- 
ties of BN nanostructures largely independent of the layer 
arrangement i 1 1 ' 13 A single layer of h-BN is thus a generic 
system for studying the optical properties of crystalline 
and nanostructured BN including the exciton-defect in- 
teraction. 

Detailed experiments on the cathodo and photo- 



luminescence of h-BN crystals and BN nanotubes have 
been performed in the recent years 2 - - — The origin 
of some of the features in the BN luminescence spec- 
tra remains unsettled. For instance, the occurrence of a 
defect-related emission band at 4 eV, composed of regu- 
larly spaced peaks compatible with phonon rcplica a 5 ' 7 ' 15 , 
can be associated with a deep level impurity. However, 
the defect that is responsible for this deep level has not 
yet been identified. Furthermore, the emission band at 
5.77 eV shows different sub-peaks on the low energy side. 
The origin of its sub-structure is unclear, since the main 
peak can be understood as being due to the recombi- 
nation of a Frenkel-type exciton 1 ^ with a degenerate re- 
combination linei 13 ' 16 Interpretations of the experimen- 
tal findings are based on different mechanisms: quasi- 
donor-acceptor pairs^, excitons coupled with structural 
defects^i 1 ^, breaking of the exact hexagonal symmetry 
(e.g., by the presence of defects) 1 ^, or a dynamical Jahn- 
Teller effect^ 7 , (symmetry breaking due to geometry re- 
laxation in the excited state). 

Theoretical works on defects in BN-structures have so 
far focused on the electronic and structural properties of 
point defects 1 ^— . The abundance (related to the for- 
mation energy) and stability of different defects has been 
addressed with detailed density-functional theory (DFT) 
calculations ! 20 ' 25 " — . The position of deep levels within 
the fundamental band gap has been calculated on the 
DFT level. However, DFT calculation suffer from the 
well-known underestimation of the band-gap. Therefore, 
besides the neglect of the strong excitonic effects in BN, 
(de)excitation spectra based on DFT energy levels may 
not be reliable enough for comparison with experimental 
spectra. 
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In the present paper we address the optical properties 
of point defects in BN going beyond the approximation 
of DFT. We focus on a substitutional carbon impurity on 
the nitrogen site, the isolated vacancies and the BN di- 
vacancy that are common defects in BN. We will demon- 
strate that defects are responsible for the emission lines 
around 4 cV. Their presence can also explain the exper- 
imentally observed splitting of the main exciton peak. 
The paper is organized as follows: in section |TT] we sum- 
marize the computation methods employed in this work; 
in section Uni we present the ground state results and the 
quasi-particle(QP) band structures; in section ITVl we in- 
vestigate the role of defects on the optical properties and 
finally in section [V] we compare our results with experi- 
ments. 



II. COMPUTATIONAL METHODS 

Density functional theory (DFT) within the local spin 
density (LSDA) or gradient corrected approximations 
(GGA) for the exchange-correlation energy has been 
widely used to study defects in semiconductors and in- 
sulators. Although DFT is in principle exact, it is lim- 
ited to the ground state and as a closed expression for 
the exchange-correlation energy functional is unknown, 
approximations like the LSDA or GGA have to be em- 
ployed. Albeit theoretically well-founded, these approxi- 
mate functionals bear inherent limitations like the well- 
known band gap problem which concomitantly affect the 
prediction of defect levels, in particular in wide band 
gap insulators. In the last years, many-body pertur- 
bation theory within the GW approximation 3 ^ has be- 
come a valid and accurate alternative to predict defects 
properties in semiconductors and insulators^— start- 
ing from the DFT ground state. Moreover, in going 
beyond the ground state, the optical properties of de- 
fects and the host insulator can be addressed by includ- 
ing the electron-hole interaction in the excited state via 
the Bethe-Salpeter equation (BSE)2±-^. 

Here we use this approach to study the role of pro- 
totypical defects in the optical properties of BN nanos- 
tructures: a substitional carbon impurity on the nitro- 
gen sublattice (Cn), the boron vacancy (Vb), the ni- 
trogen vacancy (Vn) and the boron- nitrogen di- vacancy 
(Vbn)- Calculations proceed in three steps: first we ob- 
tain the DFT ground state employing the LSDA approxi- 
mation and a plane wave pseudopotential method^ (see 
HH for details). Second we evaluate the quasi-particle 
band structure within the Go W o 30 ' 39 ' 40 approach and 
third we address the optical spectra via the BSE4i A 
very large supercell 3 ^ is employed to accommodate the 
defect and to enable at the same time the description of 
defect bound and free excitons. This becomes possible 
due to the strong localization of the bulk exciton state in 
h-BNfi2 in the layer, and tube s 13 ' 42 . 

As it has been shown that a large class of impurities, 
when introduced into h-BN, is in a spin-polarized state 



and acquires a local magnetic moment^!, we have per- 
formed spin-polarized calculations for the Cn, Vb, and 
Vn cases while this was not necessary for the di- vacancy 
Vbn that does not display a net spin polarization. The 
quasi-particle properties have been calculated starting 
from the Kohn-Sham Hamiltonian and orbitals within 
the Go Wo approximation^ 3 -: 

E n k = e n fe + Z n kReATi n k(e n k), (1) 

where AS = £ - ^ xc and V xc is the LSDA exchange- 
correlation functional, e n k are the Kohn-Sham eigenval- 
ues, and Z n k = [1 — dReAYt/dw]^ 1 is the renormalization 
factor—. The screened electron-electron interaction W 
has been evaluated within the random-phase approxima- 
tion (RPA) in terms of the dielectric function cg,G' (<7> w ) 
using a plasmon-pole model4^ Neutral excitations are 
investigated on the basis of the quasiparticle spectrum 
by solving the Bethe-Salpeter equation (see |45| for de- 
tails). Excitonic effects are analysed through a compar- 
ison of the absorption spectra for independent particles 
with absorption spectra where electron-hole attraction is 
included on the level of the BSE. 

As we are studying an isolated monolayer employ- 
ing three-dimensional periodic boundary conditions, the 
size of the quasi particle energies and the energy of the 
electron-hole attraction will depend on the inter-sheet 
distance. In Ref. [ID we had shown that the value of the 
gap of the pure h-BN sheet is 0.6 eV higher for an inter- 
sheet distance is 80 a.u. than for an inter-sheet distance 
of 20 a.u. (value which we are using in the calculations 
presented here). At the same time we had observed an 
increase of the exciton binding energy with increasing 
inter-sheet distance that almost exactly cancels the vari- 
ation of the band-gap. We expect a similar cancellation 
effect in the presence of defects but cannot exclude a 
minor effect of the inter-sheet distance on the absolute 
position of defect peaks. 

III. DFT-GROUND STATE AND 
QUASI-PARTICLE STATES 

In this section we address the electronic and geomet- 
ric structure of the four defects and then compare the 
DFT-electronic structure with the quasiparticle states 
obtained from Go Wo calculations. 

In Fig. [1] the equilibrium geometry of the Cn and the 
vacancies VnVb and Vbn is shown along with the BN 
bond distances in the vicinity of the defect center. In 
the case of subsitutional impurity all the first neighbor B 
atoms of Cn undergo a small symmetry conserving out- 
ward relaxation of 0.113 a.u. (4% of the BN-bond dis- 
tance) with respect to the ideal lattice positions. The de- 
fect retains its threefold symmetry. In the case of Vb and 
Vn the the outward displacement of the first-neighbor 
atoms is larger. More importantly, in agreement with 
Ref. [24|, we find that a (pseudo) J aim- Teller effect breaks 
the three-fold symmetry slightly. In the di-vacancy Vbn 



FIG. 1. (color online) Geometry of the defects in BN mono-layer with defects: (a) Cn', (6) Vb', (c) Vn', (d) Vbn- Defect induced 
changes of BN bond distances are indicated (in atomic units) and have to be compared with the BN distance of the perfect 
sheet (2.724a.u.). 



case, the three-fold symmetry is already reduced to a 
mirror symmetry by the defect configuration itself. We 
observed that the first N and B neighbors decrease their 
mutual distance upon relaxation, as shown in Fig. [T}i. 
The defect geometries found here are in good agreement 
with those reported in Ref. and [l8|. 

As established already in earlier work, the Cn and 
the vacancies possess localized defect levels in the band 
gap. Here we briefly review the electronic structure of 
the defects as obtained in DFT-LSDA and then compare 
it with the quasiparticle states obtained on the level of 
the Go Wo approximation where the DFT band gap-error 
is corrected. Since the prediction of defect emission and 
absorption spectra hinges on the defect level positions, a 
correct prediction of these quantities is a crucial point for 
the investigation of optical properties. The defect levels 
as obtained within the DFT-L(S)DA are reported in the 
left panel of Fig. H 

In the case of Cn the defect levels arise from the 
three bonds of the substitutional carbon impurity with 
its boron neighbors. These bonds combine into two fully 
occupied states located within the valence band and a 
third half-occupied localized state within the band gap. 
As a result of the occupation with only one electron there 
is a spin-splitting of ~ 0.5 eV between the occupied 
(spin-up) and unoccupied (spin-down) state. 

The situation is more complex for the boron vacancy. 
In this case a two-fold degenerate and a non-degenerate 
defect state are formed from the unpaired orbitals at the 
three nitrogen neighbors. In agreement with Ref. [2(| [2(| 
and [29|, these defect levels are located within the band 
gap. In the neutral vacancy the two-fold degenerate level 
is occupied by three electrons and the non-degenerate 
level remains empty. Due to the Jahn- Teller distortion a 
small splitting of the degenerate level arises in addition to 
the exchange-splitting of the half- and unoccupied states. 
Vb behaves as a triple acceptor. 

For Vn we obtain a half-filled level in the band gap, 
and two un-occupied levels almost degenerate in energy 
with the conduction bands. Due to spin-splitting their 
counter parts of opposite spin are found as resonances 
in the conduction band. Finally, in the case of the di- 



vacancy Vbn , there are four orbitals originating from the 
two boron and nitrogen neighbors. Only three are located 
in the band gap, one being occupied by two electrons and 
the upper two remaining empty. 

The above picture obtained from the Kohn-Sham 
eigenvalues gives a qualitative description of the defect 
levels in the band gap. The position of the states is, 
however, quantitatively affected by the LSDA-band gap 
error and the approximate treatment of exchange with 
this XC-functional as discussed above. Furthermore, the 
ordering of occupied and unoccupied states on the Kohn- 
Sham level may be different from that of the "true" quasi- 
particle energy levels (see e.g. H3). The Go Wo method 
yields true quasiparticle energies and an improved de- 
scription of exchange and screening effects. Due to the 
reduced dimensionality of the BN mono-layer there is an 
incomplete screening of the Coulomb interaction, such 
that the quasi-particle corrections are expected to be 
quite large. We calculated the quasi-particle band struc- 
ture for all the cases considered above. The calculated 
band gap is in agreement with the results of Ref. [H. In 
the right panel of Fig. [5] the position of the defect quasi 
particle energies within the band gap is shown schemat- 
ically. The G0W0 corrections to the LSDA levels de- 
pend on the orbital occupation and on the wave-function 
character—. The exchange splitting between occupied 
and unoccupied defect levels is strongly enhanced and the 
ordering of defect levels within the band gap is altered in 
some cases. In all cases we find that occupied levels are 
pushed closer to the valence bands and unoccupied ones 
closer to conduction bands. For Cn the spin-splitting is 
significantly increased from 0.5 eV to 3.1 eV. A large ex- 
change splitting is also present in the case of Vn- For 
Vb all un-occupied levels are pushed up by more then 
2 eV, depending on the character of the wave-function, 
and the level order is modified as well. In the case of 
Vbn, un-occupied levels are found close to the conduc- 
tion band edge, while the occupied ones remain in their 
position close to the valence band edge. This underlines 
the fact that the di-vacancy contains both donor and ac- 
ceptor levels. The general picture obtained from Go Wo 
shows the clear deficiencies of the DFT-LSDA in the dc- 
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FIG. 2. (color online) Quasiparticle energies of Cn, Vb, Vn, and Vbn as obtained from LDA and Go Wo calculations: schematic 
representation of energy position with in the band gap. The filled arrows indicate levels occupied with electrons with spin up 
or down. Dotted arrows indicate unoccupied states (holes) correspondingly. 



scription of defect level positions in the case of BN, in 
particular, for unoccupied states. Vested with the de- 
fect quasiparticle spectrum we now turn to the optical 
properties. 



IV. OPTICAL PROPERTIES 

In order to discern excitonic effects on the optical prop- 
erties, we present first in Fig. [3] the optical spectra in 
the independent particle picture (IP)^ (i.e. without the 
inclusion of the electron- hole coupling). This yields an 
overview over relevant excitation mechanisms and later 
allows one to trace down the effect of the electron-hole 
coupling. 

We start the discussion with the case of Cn- Compar- 
ing the corresponding panel of figurc[3]with the schematic 
representation of the defect levels in the band gap (Fig. [2]) 
one can identify the different absorption peaks. The first 
band at 3.4 eV is due to transitions form the extended 
states at the valence edge to the defect levels within the 
band gap; the other band at 4.4 eV involves extended 
states below the valence band edge and defect levels. 
Above 6 eV we find the bulk-like spectrum that is only 
slightly perturbed by the presence of the defect. For the 
boron vacancy Vb there are two bands at 3.5 eV and 
4.7 eV that are associated with transitions from defect 
resonances in the valence band to the localized levels in 
the band gap and a series of small peaks around 5.9 eV 
originating from valence band to defect transitions. The 
defect-related excitation bands of Vn at 4.1 eV, in con- 
trast to the previous cases, involves extended states near 
the conduction band edge as final states. Only the faint 
defect-related excitation band around 6.6 eV, not visible 
in the figure, originates from valance band states. An- 
other faint peak around 5 eV is due to transitions among 
the defect levels in the band gap. Finally, for the Vbn 
di-vacancy, we find peaks between 5.6 eV and 6.0 eV due 
to defect-defect transitions and valence-band to defect 
transitions, respectively. 



We now turn to the spectra— including the electron- 
hole interaction via the Bethc-Salpeter equation. These 
are shown in Fig. [5J The inclusion of the electron- hole in- 
teraction has essentially two effects on the defect-related 
spectra: (i) a general red-shift of the defect excitation en- 
ergies due to electron-hole attraction and (ii) a coupling 
of the independent particle excitations. This coupling 
can give rise to bound electron-hole states, the excitons. 
In our case the BSE leads to a red-shift of defect-defect 
transitions in the range from 0.6 to 1.3 eV depending on 
the character and the localization of the wavefunction. A 
common feature of all spectra is the large peak around 
5.5 eV which corresponds to the strongly localized exci- 
ton of pure h-B N 12 ' 13 ' 42 which is composed of a large 
number of independent electron-hole pairs. In the pres- 
ence of a defect, electron- hole coupling leads to defect 
bound excitons besides the "free" excitons. 




2 2.5 3 3.5 4 4.5 5 5.5 6 6.5 6.5 7 7.5 8 8.5 9 9.5 10 



eV eV 

FIG. 3. Independent-particle absorption spectra (arbitrary 
units) based on the GoWo quasi-particle band structure. In 
the left panel, the absorption cross section is increased by a 
factor 15 with respect to the right panel in order to make the 
defect-related peaks visible. . 

In the case of Cn and Vb, as one can see from Fig. [51 an 
additional small peak appears just below the main exci- 
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tons. In the case of Vbn the main exciton band has three 
sub-peaks due to the strong defect field with reduced 
symmetry. In order to disentangle the nature of the dif- 
ferent peaks, we analyze the amplitude of the electron- 
hole pairs that compose each of the exciton eigenstates A, 
as a function of the energy difference of the quasiparticle 
states of the corresponding electron and hole states: 

A\u)= \m\ 2 S(u>-E v ), (2) 

r/={eh} 

For brevity, we report A x (uj) only for Vb and Vbn (sec 
Fig. SJ. The plots for the other defects show generally a 
similar behavior. 
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FIG. 4. Amplitude of the different electron-hole pairing enter- 
ing in the BSE, for the boron vacancy and the BN divacancy, 
compared with the pure h-BN monolayer. 

Comparing for Vb the electron-hole amplitudes of the 
four excitation energies with the schematic diagram of 
Fig. [2j we can distinguish three kinds of excitations: (i) 
transitions from valence bands or defect resonances to 
the defect levels, that generate the peaks at 4.2 eV, (ii) 
the main exciton at 5.6 eV that is entirely composed of 
extended states, and (hi) a peak at 5.4 eV generated by 
the coupling of bulk excitations with transitions from 
and to defect levels in the band gap. The Cn and Vn 
cases present a structure similar to the one discussed 
above. From these results it is possible to infer three 
important consequences: (i) transition from defect levels 
to valence/conduction bands with an energy close to the 
band gap energy mixes with the bulk excitation giving 
rise to an additional peak close to the main exciton; 
(ii) transition between defect levels are renormalized by 
the electron-hole interaction but do not mix with the 
bulk excitations; (iii) the main exciton peak maintains a 
structure similar to the bulk system. 

The excitonic band of the di-vacancy Vbn is distinct 
from the other spectra as the di-vacancy possess both ac- 
ceptor and donor states of low symmetry in the vicinity 
of the band edges. The direct transitions between ac- 
ceptor and donor states give rise to the absorption peaks 
around 4.8 eV, while the resonant coupling of transitions 
between the extended conduction and valence band states 



with donor and acceptor states splits the main exciton 
in three sub-peaks respectively at 5.5 eV, 5.43 eV and 
5.37 cV (see Fig. SJ. 

The observed side peaks to the main excitonic peak 
are compatible with the fine structure observed in 
different luminescence experiments^ ' 1 5 1 1 7 . 



V. CONCLUSION 

We have investigated the role of defects on the opti- 
cal properties of BN nanostructures. We have treated 
four different plausible defects in a single sheet of hexag- 
onal BN: a subsitutional carbon impurity, a boron va- 
cancy, a nitrogen vacancy, and a BN di-vacancy. Quasi- 
particle effects strongly modify the defect-level positions 
within the band-gap and increase the exchange splitting 
between occupied and unoccupied levels. Electron-hole 
interaction (treated on the level of the the Bethe-Salpeter 
equation) not only leads to a strongly bound excitonic 
state in pure BN but also strongly renormalizes tran- 
sitions from and to defect levels. This means that the 
peaks in the optical spectra which are related to "deep 
level impurities" are strongly affected by the e-e and e- 
h interaction and can be reliably calculated only on the 
level of many-body perturbation theory (as opposed to 
the frequently used random-phase approximation using 
DFT wave- functions and energies). Experimental lumi- 
nescence spectra (as presented, e.g., in Ref. [H) display 
a prominent peak close to 4 cV with phonon replica on 
the low energy side. Comparing the four theoretical ab- 
sorption spectra (Fig. O with the experimental lumines- 
cence spectrum, the best agreement is obtained for the 
Boron vacancy where the main absorption peak of the 
low energy regime is at 4.2 eV. The Stokes shift due to 
possible relaxation of the excited state is not included in 
our calculation and limits the comparison of experimen- 
tal luminescence with theoretical absorption spectra. In 
the case of a strong (0.7 eV) Stokes shift, also the BN- 
divacancy (where a defect related absorption peak is at 
4.8 eV) could be responsible for the experimental obser- 
vation. The nitrogen vacancy (absorption peak at 4.5 
eV) might also explain the experiment. However, in this 
case a second strong luminescence peak should also be 
observable below 3 eV which seems not to be the case in 
experiments. 

The 10x10 supercell employed in our calculations 
is considerably larger than the extension of the wave- 
function of the "free" exciton in pure h-BN. Therefore, 
besides the defect-related absorption peaks, our calcu- 
lations also show the prominent absorption peak due to 
"free" excitions at 5.6 eV. This peak is doubly degenerate 
in the pure BN sheet. Coupling of the free exciton with 
one or several transitions from or to defect levels gives 
rise to a splitting of this peak and additional side peaks. 
The details of the splitting depend on the nature of the 
impurity but also on its periodic arrangement^. The 
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FIG. 5. (color online) Optical absorption spectra calculated by the Bethe-Salpeter equation . In the left panels the absorption 
cross section is increased by a factor of ten in order to make the defect-related peaks visible. Spectra are shown with a 
broadening of 0.1 eV (to simulate a typical experimental broadening) and 0.005 eV (in order to investigate the fine-structure 
of the spectrum). 



fine-structure^ of the luminescence peak at 5.7 eV (with 
additional side peaks at 5.64 eV and 5.46 eV) can possi- 
bly be explained through our calculations which clearly 
display defect-related side peaks of the "free" excitonic 
peak. Furthermore, any kind of breaking of the exact 
hexagonal symmetry can lead to a splitting of the bright 
doubly degenerate "free" excitonic peak into two peaks 
and will also (for bulk h-BN) render the "dark" exciton 
(located 0.1 eV below the bright exciton) visible^. For 
reasons of computational feasibility, we have performed 
calculations only on point defects. Other kinds of defects 
such as dislocations and grain boundaries lead similarly 
to a breaking of symmetry and/or to the formation of 
localized states and could thus as well explain the fine 
structure of luminescence spectra one would expect a 
similar interaction of the "free" exciton with dangling- 
bond states). Indeed, Jaffrcnou et al. had observed that 
the luminescence peaks at 5.64 eV and 5.46 eV occur 
mostly at grain boundaries and dislocations^ where de- 
fect states due to dangling bonds are supposed to oc- 
cur. The aim of our present study is thus not to un- 
ambiguously assign a certain type of defect to the ob- 



served luminescence spectra, but rather to show the gen- 
eral mechanism leading to the deep impurity levels and 
the formation of defect bound excitons. We also note 
that we did not study all possible point defects. A recent 
study of the chemical composition of an isolated h-BN 
sheet suggests the incorporation of oxygen atoms in the 
honeycomb lattice^. The same study also demonstrated 
that substitutional carbon atoms occur mostly in pairs, 
occupying neighboring nitrogen and boron sites. Fur- 
ther experimental and theoretical studies are necessary 
to decide if the latter two types of defects can also occur 
during the synthesis of BN nanostructures or if they are 
due to the experimental conditions in the transmission 
electron microscope where vacancies in the h-BN sheet 
are induced by the electron bombardment and may be 
subsequently filled by a substitutional atom. 

An alternative route to symmetry breaking was pro- 
posed by Watanabe et ali^l. They suggested that the 
dynamical Jahn- Teller effect (splitting of the degenerate 
exciton by polaronic effects) could explain the excitonic 
splitting in hexagonal BN. Further experimental and the- 
oretical studies are necessary to disentangle the two ef- 
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fccts. 
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